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Abstract 

We customize the simulation code FEWZ (Fully Exclusive W, Z Production) to study Z' pro- 
duction at the LHC for both = 8 TeV and yfs = 14 TeV. Using the results of our simulation for 
several standard benchmark Z' models, we derive a semi-empirical expression for the differential 
cross section, that permits the determination of Z' couplings in a model-independent manner. We 
evaluate cross sections and other observables for large classes of models, including the common Eq, 
left-right, and B — L models as a function of model parameters. We also consider a hidden sector 
Z 1 that couples to standard model fermions via kinetic and mass mixing and serves as a mediator 
of isospin- violating interactions with dark matter. We combine the results of LHC Z' searches and 
dark matter direct detection experiments with global electroweak data to obtain mass-dependent 
constraints on the model parameters. 



I. INTRODUCTION 



A simple extension of the Standard Model (SM) is the addition of an extra U(l) gauge 
symmetry, with associated neutral Z' gauge boson. Extra U(l) symmetries are a necessary 
part of many interesting new physics scenarios, including several Grand Unified Theories 
and string-inspired model constructions. Generic Z' models can have many new physics 
features, including generation-dependent couplings, Z — Z' mixing, and new fermions; see 
e.g., Refs. [THS] • We first study the simplest models with generation-independent couplings 
and no Z — Z' mixing. There are a number of such models that are theoretically relevant, 
such as the Eq GUT models and the B — L model. In Section |VI we consider a model that 



includes mass and kinetic mixing between the Z' and the Z |6j, which has applications to 
isospin- violating dark matter scattering [7]. 

We expand the simulation code FEWZ 2.1 [SJ [H] to study the production and decay of 
Z' bosons at the LHC through the process pp — > Z' — > l + l~ . FEWZ includes up to NNLO 
in perturbative QCD and is fully differential in the lepton phase space. This allows for 
the precise calculation of the Z' cross sections and differential distributions with realistic 



experimental acceptances. In Sect ions inland III , we briefly review Z' detection and introduce 



several common benchmark models that we use to demonstrate the efficacy and validity of 
our simulation. 



In Section IV we use FEWZ to derive a semi-empirical expression for the differential cross 
section , a . We show that, with sufficient data, this formula can be used to determine 

ay a cos a ' ' 

the couplings of the Z' and set limits on model parameters with good precision. 

In Section |VJ we apply our fit to the Eq model class and the production of two Z' 
bosons. The heavier of the two mass eigenstates is often assumed to be too heavy for 
collider detection, but we show that it could be accessible at the LHC for a certain range of 
mixing angles. Throughout, we focus on Z' masses of a few TeV, as the LHC lower limits 
with approximately 5 fb _1 of data fall in the vicinity of 2 — 2.5 TeV for the considered 
models 1 1 1 . 



In Section VI we study the phenomenology of a Z' scenario in which SM particles are 
uncharged under the new U(l)'. In this case, SM particles interact with a new sector 
through kinetic and mass mixing of the Z' with the Z. Such a Z' could act as a dark 
matter mediator, with isospin-violating dark matter scattering arising naturally. Then, 



there are two complementary ways to test such a model: the production of Z' resonances 
in collider experiments, and the direct detection of dark matter particles. We combine the 
data from LHC Z' searches [TTI IT2] and the XENON100 dark matter experiment [13] with 
global electroweak data to constrain the kinetic and mass mixing angles. We find that the 
electroweak, collider, and dark matter data provide comparable limits. 



II. Z' BOSONS AT THE LHC 



If the Z' couples to standard model quarks and leptons, it could be detected at the 
LHC as a resonance in the dilepton channel through the experimentally well-studied process 
pp — >■ Z' — >■ [T3HT9] . We focus on the dielectron and dilmuon channels, though decays 
to r-lepton pairs can also be useful [20]. The SM background to this process is fairly small, 
consisting mostly of Drell-Yan (DY) Z/'j* events and a smaller number of ti and multijet 
events [TT1[T2]. 

The differential cross section for the DY process is [2] 

— — = -o[Z' — > l + l~) Wz'(s, Q 2 ) + interference terms. (1) 

Cl\cZ S 

The first term on the right hand side is the pure Z' contribution, factored into two parts: a 
hadronic structure function containing the QCD dependence, Wz'(s,Q 2 ), and the partonic 
cross section, 

er(7' > 1+1-} - 1 9lL +9lR ^ 0\ 

1 j 4tt 288 (Q* - M*,y + M%,T 2 Z , ' [ ) 

The terms not shown include the interference with photons and the SM Z. For narrow 

resonances, which we define as Tz> < Q.lMz>, they are small. 

The partial width for decay into a massless fermion pair // is given by 

ri- = + (3) 

where the gf LR are the fermion couplings, which can be written in terms of the overall Z' 
coupling and the fermion charges zf. 

2 2 2 2 

9f L , R -9z>z fL R - cos2()w Z fL R . W 
We take the masses of the quarks to be negligible compared to Mz'/2. 



We assume for simplicity that the Z' decays only to SM fermions. However, in models with 
extra fermions, the Z' might also have non-SM decays. Thus, the total width is generally 
a free parameter. If the decay rate to new fermions is large, the Z' mass limits could be 
significantly relaxed due to the reduced SM branching fractions j21j . 

Dilepton resonances are not the only viable channel for Z' detection. Past work has con- 
sidered detection using decays to top quarks [22H25] and third generation fermions [26ti28] . 
weak boson pair production [29] , and weak charge measurements in atomic parity violation 
experiments [301 EE]- The latter channels are particularly useful in the case of leptophobic 
or non-universal Z' models [32]. 

Other new physics, including Randall-Sundrum gravitons or sneutrinos, could also be 
detected via a dilepton resonance similar to a Z' . There have been several discussions of 
how to differentiate such resonances from a Z' [251 E3] . 

III. SIMULATION 

We have analyzed the characteristic features of Z' production at the LHC using an 
expanded version of the simulation code FEWZ [HI E] (see Appendix). All calculations are 
done to NLO or NNLO in QCD using the MSTW2008 PDF sets [34j. Factorization and 
renormalization scales are set to \ip — — 

M z > M- 

We adopt the following standard acceptance cuts in our analysis: 

p l T > 20 GeV , 
\Vi\ < 2-5, 
\y Z r\ > 0.8 (for A FB only) . 

The first two of these cuts are the same as those applied to the muon channel by CMS and 
ATLAS [HJ H2]. The third is used to define a forward-backward asymmetry. The intial 
quark direction cannot be measured directly at proton-proton colliders. However, the boost 
direction of the Z' is preferentially in the direction of the quark, not the antiquark, especially 
for large dilepton rapidity \yu\. Valence quarks are much more likely than any other partons 
to carry a large fraction of the proton momentum. Thus, if the dilepton rapidity is large, the 
boost is preferentially in the valence quark direction. However, for small dilepton rapidities, 
the initial momenta of the quark and antiquark have similar magnitudes, so we cannot 
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use the parton distributions to distinguish between them. Therefore, placing a cut on the 
rapidity of the final dilepton system allows for a measurement of Afb |19j - 

We consider several common benchmark Z' models. The first, the sequential standard 
model (SSM), has couplings identical to those of the SM Z. It is a common benchmark 
in experimental searches. Another theoretically interesting case is the Z' B _ L , noteworthy 
because U{1)b-l satisfies anomaly cancellation conditions without the presence of exotic 
fermions or non-universal couplings [21 [351 [36] and may be a remnant of string theory [371 - 

E2|. 

A theoretically well-motivated class of Z' models derives from breaking the Eq gauge 
group via the chain [3] 

SU(5) x U{1)^ x U(l) x SU(3) e x SU(2) L x U(1) Y x U{1)^ x U(l) x . 

The new U(l) factors are associated with two neutral gauge bosons, Z^ and Z x . After 
symmetry breaking, they mix to form the mass eigenstates Z' and Z", with the mixing 
parameterized by an angle (3 [10] : 

Z' = Z x sin (3 + Z<p cos (3 , 

Z" = Z x cos (3 — Z^p sin (3 . (5) 

Since the Z' and Z" are assumed to be heavy compared to the SM Z-boson, any mixing 
with the Z is negligible and ignored here. For now, we only consider the lower-mass Z 1 for 
the three specific cases ZL for (3 = 0, Z' % for (3 = tt/2, and Z' v for tan/3 = a/3/5. We explore 
Z' and Z" detection in this class of models in Section |V] 

Finally, we consider a left-right symmetric model, Z' LR [H]. The Z' couples to the current 

J LR = a LRJm - > ( 6 ) 

where o:lr = a/ cos 2 0w9rI sin 2 0wg\ — 1 and sin 2 8\y = 0.22255 |42j. We study an example 
Z' LR with g\ = g 2 L i.e., a LR w 1.58 [H]. 

Table |l] summarizes the Z' couplings to SM fermions for these models. The overall gauge 
coupling strength is a free parameter, but it is often chosen to be consistent with a grand 
unification scenario. We follow this approach and set g\, = Cl ^2g w [18] . (This factor is not 
included in Table [l|) 

5 





z u 
L 


L 


L 


R 




4 

it 


zK 

L 


ib 

Y 


s/lO 
12 


s/lO 
12 


\/To 
12 


-VTo 
12 


— /lo 

12 


-Vw 
12 


Vw 

12 


X 


-1 
2sf§ 


-1 
2^ 


3 


1 


-3 
2^ 


1 

2^ 


3 
2V6 


7 


1 

3 


1 
3 


-1 
6 


-1 
3 


1 
6 


-1 
3 


-1 
6 


B-L 


v 


v 
6^ 


_ V J 
2^ 


V 
6^2 


V 
6^ 


_ V 

2V2 


_ V 
2^2 


LR 


1 


1 


1 


1 1 CtLJJ 


1 a LR 


1 «M 


1 


6a LR 


6a LR 




Ga LR 2 


6a LR 2 


2«Li? 2 


2a£H 



TABLE I. Z' gauge charges. For the LR model, o:lr~ 1-58 



Model 


r/M z , 


BF0U+/O 


BF(tt) 


BF(hadrons) 






0.005 


0.04 


0.12 


0.80 


0.07 


X 


0.012 


0.06 


0.03 


0.65 


0.16 


77 


0.006 


0.04 


0.16 


0.87 


0.02 


B-L 


0.014 


0.15 


0.05 


0.31 


0.23 


LR 


0.022 


0.03 


0.10 


0.90 


0.02 


SSM 


0.026 


0.03 


0.10 


0.73 


0.18 



TABLE II. Z' decay widths and branching fractions, assuming no non-SM decays. 

In Table [IT] we list for each model the total width, normalized by Mz', as well as the 
branching fractions to leptons, quarks, and neutrinos, assuming no decays to non-SM parti- 
cles. For our calculations, we use ot~ l (Mz) = 128. 

In Fig. [T] we show the shape of the dilepton mass spectrum for a variety of Z' models 
with M z , = 2.2 TeV at = 8 TeV (LHC8), and M z , = 2.5 TeV at y/s = 14 TeV (LHC14). 
The SM background is very small in the resonance mass range, so for luminosities of 20 fb _1 
for y's = 8 TeV or 100 fb _1 for x/s = 14 TeV, these resonance peaks should be clearly 
distinguishable. Distributions are calculated at next-to-next-to-leading order. We neglect 

6 




FIG. 1. Simulated NNLO dimuon invariant mass spectrum for a variety of Z' models with Mz> = 
2.2 TeV and y/s = 8 TeV (left) and M z > = 2.5 TeV and y/s = 14 TeV (right), both with a 
luminosity of C = 100 fb _1 . 

detector effects such as energy resolution smearing, which can be significant for narrow Z' 
widths (T/Mz' < 0.1). To extract the mass and width of the Z' from experimental data, 
these effects would need to be taken into account by convolving with a Gaussian resolution 
function. For ATLAS, the resolution width is about 1% for electrons and 5% for muons [TT] . 
For CMS, the resolutions widths are 1 — 2% for electrons and 4 — 7% for muons [12] 

By integrating over the peak region, we can determine the cross section as a function of 
Mz> for each model. In Fig. |2j we show the mass dependence of the integrated peak (±3IV) 
cross section for our six model examples at NLO. Note that the mass dependence is largely 
model-independent. 

IV. ANALYSIS 

We now present a few empirical formulas from our simulations, which allow the cross sec- 
tion and the differential distributions over a wide range of model parameters to be predicted. 
They are also useful in extracting coupling information from experimental data. 

The LHC production cross section of the Z' depends on the mass, width, and couplings. 
This dependence can be parameterized in the narrow width approximation by 

a(pp ^Z'^ /+n) = [p (<4 + gl R ) + (1 - p) (<4 + «4)] B(Z> Z+r) f(r z >) . (7) 
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FIG. 2. NLO cross section for the process Z' — > l + l versus Z' mass, integrated over the dilepton 
invariant mass peak region (±3r) for yfs = 8 TeV (left) and -y/s = 14 TeV (right). 

The parameter p is model-independent. It quantifies the fractional contribution to the cross 
section from up-type quark events. The dependence on the Z' mass is contained in f(r Z i). 
Following Ref. |3], we use the empirical representation, 

i Wi 1 " 



f irz ,) = M--i) l-l • («> 



of the structure function dependence, where r z > = 



As in Refs. [21 [18], we define the quantities 



M' 

6q = 24^fc ( ^ _ & )( & _ ^ • (10) 
In this notation, the cross section formula takes the simple form: 

a(pp -> Z' -> / + r) = [pc u + (1 - p) c d ) f(r z ,) . (11) 

We use a NLO FEWZ simulation of the six model examples in fits of the parameters of 
Eqs. Q and Q. First, we fix Tz> and fit the simulated cross sections for the six models to 
the form Ac u + Bed, then take p = -A^- We repeat this procedure for several values of r z > 
over the range, 0.1 - 0.6 (M z > = 1.0 - 4.5 TeV for = 8 TeV; M z > = 2.0 - 7.0 TeV for 
a/s = 14 TeV), to determine how p changes with r z >. We find that this dependence can be 
approximated by the function 

p(r Z /) = 0.77 - 0.17tan" 1 (2.6 - 9.5 r z >) ■ (12) 



p 

1.0 r 
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FIG. 3. The parameter p (as defined in Eq. j7j) as a function of rz 1 = ~§~> along with an approximate 
fit to an arctangent function. Notice that simulated values (points) are shown for both 8 TeV 
(squares) and 14 TeV (circles) and agree well. 
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FIG. 4. FEWZ simulation (points) and fitted curve (lines) of the integrated Z' cross section for 
y/s = 8 TeV (left) and y/s = 14 TeV (right). We show only the Eq tj, left-right, and sequential 
models for simplicity. 

The result is shown in Fig. [3] 

To determine f(rz>), we fit the normalization A + B to the form given in Eq. ([8]). We 
find that it is possible to quantify the dependence on the PDF's within 20% over the entire 
range of masses accessible at the LHC and all six trial models with 

17.5 / „ \ -15.0 

I 3200 fb ( ^- - 1 ] 

f{rz> 



3200 fb 
43.3 fb 



l 

r z' 

1 

r z' 



16.8 



1 

r z' 



-13.1 



y/s = 8 TeV , 
l| ( ,'- ) yfi=UTeV. 

A comparison of our fit to the simulation is plotted in Fig. |4} 

After a Z' is detected, the next step will be to measure its couplings. Measuring couplings 
with precision requires a large number of events, which may not be the situation for a very 
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heavy or weakly coupled Z' . However, with enough events, an accurate measurement of c q 
and e q is possible, as we demonstrate by considering a set of simulated measurements of a 
2.5 TeV Z' &ty/s = 14 TeV. 

The coupling combinations c q and e q of Eq. ^ can be determined by considering the 
differential cross section, integrated over the resonance peak: 
da in , r M z> +3F da 



dy d cos 9 Jm ,-3T dy d cos 9 dQ 



dQ 



3/8 (1 + cos 2 9) [p c u h\{y) + (1 - p) c d hf(y)] f(M z ,) (13) 
+ cos 9 [pe u h u 2 (y) + (1 - p) e d h d 2 {y)} g(M z ,) . 



The functions h\ 2 are normalized so that integrating over cos 9 and y yields Eq. ( 11 ). g(Mz') 
represents the mass dependence of the cos# term, similar to f(Mz>). For Mz> = 2.5 TeV, 
its value is g(Mz>) = 506 fb. For the same mass, we find f(Mz') = 1050 fb. 

We have already determined p and f(Mz>), so all that remains is to find h\ 2 . In order 
to separate the up and down contributions to the differential cross section, we define two 
distinct scenarios in which the Z' couples exclusively to one type of quark (and to leptons). 
Then, for each scenario we simulate two differential distributions. The first is the dilepton 
rapidity distribution, which allows us to determine h\. To determine h q 2l instead of using 
the total cross section, we consider the quantity d ^ F d ~ B ^ , where F is the number of lepton 
pairs scattered in the forward (cos 9z> > 0) direction in the Collins-Soper frame jl3] and B 
is the number scattered in the backward direction. The resulting distributions are easily 
distinguishable, as can be seen in Fig. [5} In addition to the simulated NLO data, Fig. [5] 
shows approximate curves h\ 2 (y) . The specific curves used to approximate the normalized 
distributions at 14 TeV are 

0.59 



K(y) 



e 4.6(j/-0.84) + I ' 



K{y) 



K{y) 



e4.6fe-0.63) _|_ I ' 

2.5 (1 - e-°- 60f ) 

e 5.7(j/-0.84) + I 

11 (1 - e" - 20 ^ 



h i(y) = 75.8(Vo.63) + ( ■ 

The rapidity distribution for an arbitrary Z' is a linear combination of h u and h d - 



By fitting observed data to Eq. (13), one can determine the four coefficients c q and e q . To 



demonstrate the feasibility of this method and estimate the error, we treat our simulation as 
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FIG. 5. Normalized differential cross section versus the dilepton rapidity yu for Mz 1 = 2.5 TeV at 
LHC14. Points are simulated data and lines are fitted curves. 
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FIG. 6. 95% C.L. regions for the couplings c u and for an average experiment (statistical errors 
only) for C = 100 fb -1 (left) and C = 1 ab _1 (right). The values of c u and Cd for the family lie 
on the dashed red contour. Points are the theoretical values. 



a pseudo-experiment and extract c g and e q for the models discussed in Section III Figures [6] 



and [7] show the 95% confidence level (C.L.) regions for our example case of a 2.5 TeV Z' 
at 14 TeV. We see that 100 fb^ 1 allows for some model differentiation, while for 1 ab _1 of 
luminosity the confidence regions are narrow. In the right panel of Fig. |6j the red dashed 
contour shows the values of c u and q for the E e models as a function of the mixing angle j3, 
from which we can see that some model differentiation should be possible. The tilt of the 
ellipses arises from the requirement that the up and down quark contributions are summed 
to give the total cross section, restricting c u and Cd to lie on a line. 
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FIG. 7. 95% confidence level regions for the couplings e u and e d for an average experiment for 
C = 100 fbr 1 (left) and C = 1 atr 1 (right). Points are the theoretical values. 

Within the Eq model class, our fit can also be used to place limits on j3 using the least 
squares method. Using Eqs. ^ and rt9J), we can write Eq. (13) in terms of (3 instead of c q 
and e q . At the 1-sigma level, we find 



0.44vr < /3 X < 0.54tt , 

0.97vr <fy< 1.03tt or 1.14vr < < 1.17tt 
0.20tt < (3 V < 0.22tt . 



(15) 



If the number of events is low, an analysis can still be done by integrating over the 
distributions h\ 2 , leading to a system of four equations in four unknowns that can be 
inverted ITS]. Note that the results in Ref. US] account for both statistical and PDF errors, 

E 1 E 3 ' 

while ours only include statistical errors. 

From our fit, we can also determine the forward-backward asymmetry, 

F-B 



A 



FB 



F + B 



(16) 



On the Z' peak, Afb depends solely on the Z' couplings to fermions. Using Eq. (13) 
integrated over appropriate ranges of cos 9, we can write 

(p«2 e u + (l-p) a d 2 e d )g(M z >) 



A 



FB 



(pa u lCu + (l-p) a d lCd )f{M z , 



(17) 
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FIG. 8. Predicted on-peak forward-backward asymmetry versus mixing angle (red curve) and the 
error bands for the FEWZ simulated values for Mz' = 2.5 TeV and \fs = 14 TeV for the Eq models 
iff, X; and r] parameterized by angle /3 (left) and the LR, B — L, and \ models parameterized by 



the angle a; see Eq. (18) (right). 



where a? 2 = 2 f ymax h^ 2 (y)dy. We choose y ma x = 2.5 and y m i n = 0.8, the rapidity cuts 
discussed in Section IIII1 

In the left panel of Fig. [8j we show the simulated uncertainties on Afb for three Eq 
models with 100 fb _1 of data. The red curve shows the predicted values for Afb versus 
the mixing angle (3 defined in Eq. (pj). The error bands can be used to set limits on the 
value of (3, though for < (3 < 0.37T, the Afb values are all very small and difficult to 
distinguish. The couplings of the LR, B — L, and x models can also be parameterized (up 
to a normalization factor) by an angle a with [Tj 



Z' = cos aZ-y + sin a Zy . 



(18) 



This is equivalent to 



Z' LR = cos 9 LR (-Z B ^ L ) + sin LR Z R 



(19) 



where arctan(a £fl ) = 8 LR = a + arctan ^/2/3. As for the three E e models, we show A F b 
versus a in the right panel of Fig [8] 

Off peak, the Z' interference with the SM background contributes to the asymmetry, so 
Afb varies significantly with dilepton mass. Figure [9] shows the effect of this interference 
on Afb- Both the shape of the curve and the peak value are highly model-dependent. 
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FIG. 9. Simulated forward-backward asymmetry for a number of Z' models with Mz> = 2.5 TeV 
and y/s = U TeV. 

V. E 6 MODELS 



We now consider Eq grand unification scenarios in more detail. An Eq gauge group can 
be broken down into either a rank-5 or rank-6 subgroup. In the rank-5 case, this leads to one 
additional Z' boson, the Z' discussed above. In the rank-6 case, there are two additional Z's, 
corresponding to the additional U(l)$ and U(l) x groups in Eq. (pjj). The mass eigenstates 
are Z' and Z" of Eq. @. We justifiably ignore small mixings of the Z' bosons with the 
SM Z. Often, the Z" is assumed to be very heavy, leading to an effective rank-5 group, as 



was the case in the models we considered in Section III In Fig. 10, we show the branching 
fractions and total width of each additional boson as a function of the mixing angle (3. 




FIG. 10. Total width, normalized by the mass, and branching fractions to lepton, top, and neutrino 
pairs for the lower mass eigenstate Z' (left) and the higher mass eigenstate Z" (right). For both, 
we include only decays into SM fermions. 
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FIG. 11. Integrated peak cross section for the process pp — > Z 1 — > I' 
angle /3, with M z > = 2.5 TeV at LHC14. 

In the rank-6 case, the masses of the Z' and Z" are related by 



as a function of mixing 



cos (3 + v 15 sin (3 \ /cos/3 



(20) 



M Z " J \ V15 cos (3 - sin (3 J \ sin (3 

This relation assumes that the U(l)' symmetry breaking scale is much higher than the 
electroweak scale. Since the experimental lower bound on the Z' mass is currently about 
2 TeV for Eg models, and there are tight limits on mixing with the SM Z, this is a justified 
assumption. 

Requiring that the Z" be heavier than the Z', and that both masses are positive leads to 
the condition 

- Vl5/A < cos/3 < (21) 



Notice that the Z' with cos /3 = a/5/8 is excluded from the range in Eq. (21 ). Additionally, 
both the ZL and Z' , with cos/3 = 1 and cos/3 = 0, respectively, have M z > <C M Z u . Therefore, 
we would not expect the LHC to detect a heavier mass eigenstate for the three E s models 
considered in Section I 



Using the empirical equations determined in Section IV, we can now calculate the cross 
section for the Z' and Z" . Both cross sections are a function of just two free parameters: 



M z and (3. In Fig. 11, we show the integrated peak cross section for Z' production and 



decay at LHC14 versus mixing angle. In Fig. 12, we plot the same quantity for the Z" over 



the allowed range of mixing angles. Here we see that if a Z' were to be discovered in the 
mass range of 1 — 2.5 TeV or so, the higher mass Z" could be accessible at the LHC as well, 
within a certain range of mixing angles. 
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FIG. 12. Integrated peak cross section for the process pp —> Z" — > at LHC14 versus mixing 
angle (3. The mass of the heavier eigenstate is determined by the Z' mass and the mixing angle. 

VI. DARK MATTER INTERACTIONS 

There are a number of scenarios in which a Z' boson can serve as a dark matter media- 
tor [13H1H]- These models often include mixing between the gauge bosons, leading to small 
effective couplings between dark matter and SM fermions [HI • We study this possibil- 

ity, paying particular attention to the possibility of isospin-violating dark matter scattering, 
which occurs naturally in the case of a Z' mediator. 

We start by considering the specific case of a single new U(l)' symmetry, under which 
the dark matter is charged, but the SM particles are not. The corresponding Z' boson mixes 
with the SM Z through both kinetic and mass mixing terms. The Lagrangian in this case 
is [5U] 

£ =£>sm — -^Z'uvZ' + -MpZ 1 fj,Z' 

- 9' E ^"Ov ~ ~ ^Z'^ + 5M 2 Z>,Z» . 

i 

Here sine and SM 2 parameterize the kinetic and mass mixing between the Z' and the Z. 
As usual, B^, W^y and Z' ^ are the field strength tensors for U(1)y, SU(2)l and U(l)'. ipi 
are the fermion fields, and fy and f\ are the vector and axial charges of the fermions under 
U(l)'. For simplicity, we consider the case where all SM fermions have f v — f\ — 0. This 
choice leads to rather weak couplings between SM fermions and the new Z 1 ', which avoids 
current LHC bounds on Z' production, fy must be nonzero to allow for spin-independent 
scattering of dark matter on nuclei. 
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We define two additional parameters for convenience: 



2 



2 ' 



tan 2f - -2 cos + 5^ sine) 

5 M|,/M| - cos 2 e + sin 2 e + 26 sin e ' 1 J 

The Z' couplings to DM and the shifts in the Z couplings are proportional to £. Since £ is 
approximately proportional to M7,? for 5, e <C -g^-, these couplings must be small for heavy 
Z'b. 

The physical states A^, Z^, and Z'^ are obtained through two sequential transformations. 
First, we diagonalize the field strength tensors; then, after SU(2) x U(l) breaking, we 
diagonalize the resulting mass matrices. After these transformations, the physical states are 
related to the original (unhatted) states by 

A ll = A li + c w sin e Z'^ , 

Z^ = cos £ (Z'n — s\y sin eZ') + sin £ cos eZ'^ , (24) 
= sin e cos — sin £,(Z^ — sin eZ^) . 

We can write the couplings of the physical states to fermions in terms of the oblique 
parameters, S, T, and U [50] : 
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(25) 



where 



and 



sl = s 2 w + 2 - 2 (jaS - (? w s 2 w aT) , (26) 

s 3 1 1 

■s = % + 2 _^ 2 (v^-aS 1 - -aT) . (27) 
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FIG. 13. Left: The 90% C.L. allowed region for S and T based on the global data for Mh = 
117 GeV [12]. Right: This region translated into a 90% CL region for e and 5 for M z > = 1000 GeV. 



The contributions to S and T due to the Z' are, to second order in £ [50] . 

aS = A^c^ySw tan e , 

«T = e 2 (^-l) + 2£ SH/ tane. 

S and T are constrained by fits to the global electroweak data, as shown in Fig. [13 
The best fit values for an assumed Higgs boson mass of M H = 117 GeV are 



(28) 
(29) 



S = 0.03 ±0.09, 
T = 0.07 ±0.08, 

with a strong correlation of 87%. The S and T values would change very little for a Higgs 
mass of 125 GeV, as may be suggested by recent LHC observations [SSI EE]- The constraints 
on S and T can be translated into constraints on the Z' mixing angles (and consequently 



the couplings and cross sections) with a simple Monte Carlo. The result is shown in Fig. 13 



Note that the kinetic mixing angle e can be quite large, assuming that 5 is small enough. 

By sampling values in the region of allowed e and 5, we can determine the range of possible 
values for the dark matter scattering cross section and the LHC Z' production cross section. 
These can then be compared to the results of experimental searches to further restrict the 
mixing angles, Z' mass, and f n / f P (defined below). 
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FIG. 14. Contour plots of f n / f p and f p /f n as a function of e and 5. No limits are placed on the 
values of e and 5. 



For direct detection experiments, we are interested in interactions between dark matter 
and atomic nuclei [6]. To determine this cross section, we start with an effective dark 
matter-quark coupling given by 



V,A V,A V,A V,A 

h v,A = 9 x z'9 fZ ' 9 x z 9 fZ 
f ' Ml, + M 2 7 



which leads to DM-nucleon couplings of 



fn = 2b V d + bl , 

f P = b v d +2b v u . 



(30) 



(31) 



Using these expressions, we can determine f n f f p as a function of e, 5, and Mz>- Contour 
plots of f n / f P and f p / f n as functions of e and 5 are shown in Fig. 
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Comparing with 



Fig. 13, we see that S and T place no limit on the value of f n /f p , though the mixing angles 



are more tightly constrained for some values of f n jf p than others. Therefore, the limit on 
the dark matter scattering cross section will vary significantly as a function of f n j f p - 

Finally, we can write the dark matter-nucleus spin-independent scattering cross section: 



V P V fp 
where a p is the spin-independent DM-proton cross section, 



p ■ 



(32) 



0V. 



64tt 



(33) 
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FIG. 15. The 90% C.L. upper bound on the spin independent dark matter-proton scattering cross 
section from XENON100 as a function of f n / f p . 

and Z and A are the atomic and mass numbers of the detector material, is the reduced 
mass of the dark matter-nucleus system, and \i p is the reduced mass of the dark matter- 
proton system [7J. 

It is common to present the spin-independent cross section for a nucleus with Z protons 
in terms of the cross section for scattering off a single nucleon, making the assumption that 
fn — fp- To compare a model with data, we must also account for possible isospin violation. 
To convert the nuclear cross section to a proton cross section a p , we multiply by [7J 



Y Ji w 2 A \z + {A l -z)f n /f P Y 



(34) 



This factor is derived from Eq. (32) by summing over all stable isotopes of atomic number Z, 



weighted by their natural abundances r/i. In Fig. 15, we show the effect of isospin violation 
on the xenon (Z = 54) cross section, normalized to the current limit from XENON100 
a p < 10~ 44 cm 2 at f n /fp — 1 and M x = 100 GeV [T3|. Isospin violation can relax the 
bound by several orders of magnitude, with the least stringent bound occuring around 

U/f p = -0.7 ezj. 

With the XENON bound generalized to all values of f n J f p , we can use it to place limits 
on the model parameters. We start by choosing a random sample of points (e, £) within 



the allowed region shown in Fig. 13 and then calculate the proton SI cross section for each 



point. In Figs. 16 and 17, we show the dependence of the cross section on f n /f p and M z >, 



respectively. From Fig. 17, we see that the largest cross sections occur for f n /f p ~ 0.35, 



while the lowest occur for /„/ f p < —0.5, which is where the XENON bound is the most 
relaxed. For both figures, we have set g x — 1 and M x = 100 GeV, which in general are free 
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FIG. 16. The cross section fro DM scattering on protons, a p versus Z' mass for f n /fp = 1- 
The thick red line at o v = 10 -44 cm 2 is the XENON100 limit [13J; dots are the cross sections 
corresponding to pairs (e, 5) sampled uniformly over the allowed values. The dark matter mass 
and coupling to the Z' are set to M x = 100 GeV and g x = 1, respectively. 




FIG. 17. o v versus f n /f p mass for M' z = 1000 GeV. The thick red line is the XENON100 
upper bound; dots are the cross sections corresponding to pairs (e, 5) sampled uniformly through 



the allowed region shown in Fig 13 The dark matter mass and coupling to the Z' are set to 



100 GeV and g x = 1, respectively. 



parameters. However, they have no impact on the qualitative features of the distributions, 
as they only enter the overall normalization factor. They influence our ability to place limits 
on the other parameters; in particular, the g x dependence of the cross section means that 
all constraints can be evaded by choosing a small enough coupling. For Mz> = 1000 GeV 
and M x = 100 GeV, the XENON bounds are evaded for all values of f n /fp with g x = 0.58. 

Collider searches can also help constrain the mixing angles for a dark Z', as the Z' 
acquires small couplings to SM fermions via mixing effects. Z' production and decay to 
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FIG. 18. a(pp — > Z' — > l + l ) versus M' z . The thick red line is the most recent ATLAS limit 
dots are the cross sections corresponding to pairs (e, d) sampled uniformly through the allowed 



region shown in Fig 13 



leptons at the LHC depends primarily on the mixing angles and the Z' mass, with only a 
small dependence on the dark matter properties through the decay width of the Z'. Using 



the couplings in Eq. (25), we can apply the analysis of Section IV to calculate the cross 



section for Z' production at the LHC. We find the ATLAS predictions of the cross section 
for the various Z' models given in Fig. 2 of Ref. [11] are well parameterized by 



<j B = 2200 fb ( 1 

Jz 
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1 

r z < 



-12 



\pc u + (1 -p) c d ) 



(35) 



With this equation, we can use the current ATLAS limits to restrict the parameter space. 
Since there are no direct couplings between the Z' and the SM, the lower bound on the Z' 



mass is much less stringent than for the models considered earlier. In Fig. 18, we show the 
upper bound on a(pp — > Z' — > set by the 5* and T parameters for a/s = 7 TeV. 

Finally, we can combine the limits from XENON100 and ATLAS to constrain the model 



parameter space in terms of e, 5, and Mz>- The results are shown in Fig. 19 , again with M- 



100 GeV and g x — 1. The electroweak, dark matter, and LHC data provide complimentary 
bounds, with S and T more strongly limiting the degree of mass mixing, while XENON and 
ATLAS provide more stringent bounds on kinetic mixing. The bounds on 5 and e relax as 
M z > increases and the Z and Z' decouple. 

The invisible decay width of the Z and the muon anomalous magnetic moment have also 
been used to constrain the mixing parameters e and S [51] , but they are less restrictive than 
S and T for M z > in the range 300 - 1500 GeV. 

If Mz > M x , the invisible decay width of the Z must be considered. The decay width is 
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FIG. 19. Allowed region for e and 5 for M z > = 500 GeV (left) and M z > = 1000 GeV (right). 
The dotted blue lines represent the limit from global electroweak data, the solid red lines are the 
limit from XENON100, and the dashed green lines are the limits from ATLAS at y/s = 7 TeV for 

£« l flr 1 . 



proportional to £ 2 , so the experimental 1.5 MeV bound [42] is avoided as long as £ < 0.95. 
There are also corrections to the muon anomalous magnetic moment [5T] . 



8a 



/! 3irc 2 w s 2 w M 2 Z , • 



(36) 



However, the contributions from the Z' are very small, because of the dependence on 
and £ 2 . The current experimental limit is Sa^ ~ 3 x 10~ 9 jl2], which is several orders of 
magnitude larger than the expected contribution except for very small Mz>- 



VII. SUMMARY 

Our customization of the FEWZ simulation code allows for extensive studies of Z' pro- 
duction and decay at the LHC at NLO and NNLO. Using the results of our simulation for 
representative benchmark models, we derived an empirical formula for the double differential 
cross section , d \ a Q . This formula can be used to study broad classes of models easily and 

ay a cos j j 

to determine the Z' couplings to fermions without prior knowledge of the underlying model. 
In a model-dependent analysis, it can also be used to set limits on model parameters. In 
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the case of the i? 6 -derived models, without accounting for systematic uncertainties, we find 
that the mixing angle /3 may be determined within 0.17T with 1 ab _1 of data for a mass of 
2.5 TeV. For E 6 models we also showed that within a range of mixing angles, two extra 
neutral gauge bosons should be within the reach of the LHC. 

Finally, we considered a more general Z' model with kinetic and mass mixing, which has 
interesting implications for dark matter detection and hidden sector theories. Even if the 
SM is uncharged under a new hidden sector U(l)', mixing could induce couplings strong 
enough that the Z' could be produced at the LHC and mediate dark matter scattering on 
nuclei, without violating limits from global electroweak data. In the case of a heavy Z' 
[Mz> 3> Mz) we find that the limits on the S and T parameters can be combined with 
XENON100 and LHC data to restrict the range of allowed mixing angles. For a relatively 
light Z' at 500 GeV, the mass and kinetic mixing parameters S and e must both be less than 
about 0.2. For heavy Z's, these mixing parameters are unrestricted. 
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Appendix: FEWZ Customization 

The FEWZ code allows for an extensive analysis of the SM Drell-Yan process, including 
NNLO effects and the influence of phase space cuts. With a few modifications, its features 
can be used to study Z' production as well. The details of the original simulation are 
provided in Refs. [8j |9]. 

The QCD factorization theorem allows us to write the Z' production cross section in 
terms of the partonic cross section and the proton structure functions as follows: 



Since the addition of a Z' affects only the partonic cross section, and not the structure 
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functions, we do not need to alter the Monte Carlo portion of FEWZ. 

The parameters of the Z' model are specified in the input file for each run. For each 
model, the user sets the mass, total width, partial width to leptons, and couplings of the 
Z' . The input file also includes the SM parameters and kinematic cuts. We also include a 
switch to turn on and off the Z' contribution, so that calculations of the SM background 
can still be done. 

The Z' parameters are read into FEWZ and used to calculate "weights" (related to the 
partonic cross section of various subprocesses) for the integration routine. Preserving the 
structure of the SM calculation, we insert additional Z 1 contributions to the partonic cross 
section. This includes the interference terms mentioned in Eq. ([I]): 

v , _ 9z>9x ( zi L z? L + z lR z* \ (Q 2 - M 2 ZI ) (Q 2 - Af 2 ) + M z , M x Y z , Y x 

(Q*-M x ) 2 + M x r x ~\ ' 

where X — j, Z. Hereafter, the Monte Carlo integration proceeds without alteration. 
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